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ABSTRACT 



^ ■ Production of rare particles within rapidity gaps has been proposed as a back- 

ground-free signal for the detection of new physics at hadron colliders. No 
complete formalism accounts for such processes yet. We study a simple lowest- 
order QCD model for their description. Concentrating on Higgs production, 
we show that the calculation of the cross section pp — > ppH can be embedded 
into existing models which successfully account for diffractive data. We extend 
those models to take into account single and double diffractive cross sections 
pp — > HXiX 2 with a gap between the fragments X± and X 2 . Using conserva- 
tive scenarios, we evaluate the uncertainties in our calculation, and study the 
dependence of the cross section on the gap width. We predict that Higgs pro- 
duction within a gap of 4 units of rapidity is about 0.3 pb for a 100 GeV Higgs 
at the Tevatron, and almost 2 pb for a 400 GeV Higgs within a gap of 6 units 
at the LHC with y/s=U TeV. 
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1. Rapidity gaps as a detection tool 



At hadron colliders, rare high-Q 2 processes can be produced in the usual hard collision 
picture, where two partons, one from each nucleon, collide head-on and annihilate into the 
rare particle one is trying to produce. The big drawback of such a production mechanism 
is that hadrons cannot simply loose their partons, since the annihilation of two partons 
leaves the initial nucleons in a coloured state. This means that they will not fragment 
independently, but rather that a colour string will link them. When they pull apart, 
the spectator partons will break the string. Hadronic matter will then populate most of 
the central phase space. However in roughly 10% of the hard scattering events at HERA, 
further interactions neutralize the colour of the proton fragments, and lead to large rapidity 
gaps [§. 

If one manages to produce a rare particle without changing the colour of the nucleons, 
then very little hadronic activity should be present in the event. Of course, both nucleons 
will still fragment, as the production of a heavy particle imbalances the kinematics of their 
partons, but the two fragmentations will be independent, as they are not correlated by 
colour strings. Hence the produced hadrons will follow the direction of the initial nucleons 
and there will be a large rapidity gap with no hadronic remnants. The rare particle will 
often be produced within the gap, and hence its decay should be essentially background- free. 
Bjorken [§] has suggested using such rapidity gaps as a means of detecting new physics. 

Colour-singlet exchange between protons is in fact a very common event, and it has been 
known and observed for a long time. Indeed, the simplest kind of event with a rapidity gap 
is the elastic scattering pp — > pp, which accounts for about 20% of the total cross section at 
the Tevatron || [|. Another 20% of the cross section is due to single diffractive scattering 
pp — > pX, and to double diffractive scattering pp — > XiX 2 , both processes leading to a 
large rapidity gap between produced hadrons. Hence the study of rare particle production 
in a rapidity gap can be seen as the inclusion of a high-mass component into soft, low-Q 2 
diffractive physics. 

Several attempts have been made to describe this process within a structure function 
formalism, following Ingelman and Schlein ||, and the demonstration by UA8 that hard 
diffractive scattering does exist ||. These works use a pomeron structure function, and 
treat the pomeron in a way similar to a photon. The major drawback of this approach is 
that we do not know what the pomeron is made of, and whether the concept of structure 
function holds for it. Similarly, the pomeron structure function is not well measured, even 
if it exists, and its flux factor is also subject to question 0. 

The cleanest estimate of Higgs production is presumably that due to Bialas and Land- 
shoff [|Tt| , where they perform a two-gluon exchange calculation. However, their calculation 
suffers from a few drawbacks: 

• They rely on the Landshoff-Nachtmann ]T^] model of the gluon propagator, which is 
taken to be a falling exponential. Although this model has reasonable phenomenological 
support, it is not adapted to higher-order calculations; 
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• In the case of diffractive scattering, i.e. when the proton breaks up, they rely heavily on 
Regge theory, as their cross section converges only because of the non-zero Regge slope; 

• They totally neglect exchanges involving several quarks in each proton ; 

• They neglect the effects of the longitudinal parton kinematics, which is essential for heavy 
particle production. 

In this paper we remedy all the above problems: we introduce new form factors, which 
describe diffractive scattering, and explain how the longitudinal kinematics comes into play. 
We keep the exact kinematics of the problem, and show that the cross section are IR finite 
even in perturbative QCD. We then estimate the size of the cross section, and obtain a 
result surprisingly close to that of Bialas and Landshoff. 

It is well known that lowest-order QCD produces surprisingly good estimates of diffrac- 
tive processes. In the present case, we tune the form factors so that they reproduce the 
measured elastic and diffractive cross sections. We then embed Higgs production within 
soft cross sections. 

Our first task will then be (in Section 2) to review the calculation of the elastic cross 
sections || |j , and to extend it to single- and double- diffractive processes. We will then (in 
Section 3) lay out the formalism which embeds high-mass particle production in diffractive 
physics. Although the method is quite general and can be used for a variety of processes, 
we shall concentrate in Section 4 on Higgs production. We show that our treatment leads 
to a cross section which is about 20% higher than the estimate of Bialas and Landshoff, and 
that the details of the infrared region do not play an overwhelming role in the determination 
of the total rate. We also find that the evaluation of the elastic cross section depends very 
little on the details of the model. We then evaluate the inelastic diffractive contribution 
pp — > XiX 2 H, which turns out to be a factor ~ 3 larger than the elastic one. The final 
section contains a summary of our results and conclusions. 



2. QCD models of diffractive physics 

Colour-singlet interactions between nucleons is an old topic. Regge theory attributes these 
interactions to the exchange of mesons, grouped into Regge trajectories according to their 
quantum numbers. In the high-s limit, one expects the hadronic amplitude to be a sum 
of simple poles, each pole corresponding to the exchange of the particles lying on a given 
Regge trajectory: 

i(M) = EA 2 (^ a, %»W) (2-1) 

i 

with £ the Regge signature factor. A(s, t) is normalised so that the elastic cross section is 
given by da e i/dt = \A\ 2 /16tt 2 , which through the use of the optical theorem gives a total 
cross section a to t = Im (A(s, t = 0)). 

The leading meson trajectories are the (degenerate) trajectories of the p and / mesons, 
clearly present when plotting the meson states in a J vs M 2 plane. When continued to 
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Figure 1: The lowest order Feynman diagrams for pp colour-singlet exchange 



negative values of M 2 , these trajectories are responsible for the fall-off of the total cross 
section at small s, a to t ~ 1/ \fs. At higher energies, the cross sections grow, and the most 
natural assumption is that another term in Eq. ( |2.1| ) is responsible for that rise although 
there is no observed trajectory to guide us. Hence one new Regge trajectory was invented, 
the pomeron, which is whatever makes the total cross section rise at high energy (at 1800 
GeV, more than 99% of the total cross section is due to pomeron exchange). In modern 
day QCD, non-mesonic exchanges are attributed to multi-gluon exchanges, which naturally 
produce a rising elastic cross section in perturbation theory. 



Several other properties of high-s colour-singlet exchanges can be inferred from data [IT]: 

• The pomeron trajectory of Eq. ( |2.1| ) is well fitted from ^ = 5 GeV to 1800 GeV, and 
from t = to — 1 GeV 2 to the form 

a P (t) = 1.08 + 0.25* (2.2) 

• The pomeron coupling to protons, /3p(t), is well approximated by the coupling of photons, 
i. e. by the Dirac elastic form factor: 

F (t) - ( 3 ' 53 - 2 ' 79t ) ( 2 31 

ll } ~ (3.53 -f) (l-t/0.71) 2 1 - 6) 

• Colour-singlet exchanges are factorizable: the ratio of the differential cross section for 
pA — > XA to that for pA — > pA does not depend on A [Q. Another related property is 
the quark counting rule: hadronic cross sections are proportional to the number of valence 
quarks contained in the hadron. This is tested not only in irp vs pp cross sections, but also 
in cross sections involving strange quarks, e.g. 2a tot (Q~p) — <j tot {T l 'p) can be predicted 
from a tot (pp), a tot (Tcp) and a tot (Kp) fL3f. 



2.1. QCD models for elastic and total cross sections 



The simplest QCD colour-singlet exchange consists of two gluons [T^]. and is shown in 
Fig. |l|. If we exchange these gluons between quarks, we obtain the following (IR divergent) 
amplitude B: 



A q (s,i) 
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(2.4) 
(2.5) 



Here bold face variables represent transverse momenta. The transverse momenta k a and 
k& are the components of the gluon momenta transverse to the direction of the quarks and 
A is the total momentum exchanged by the quarks (A 2 = t). We have introduced two 
gluon squared masses a a and a&. 

The quark counting rule suggests that we describe a proton as made of its 3 valence 
quarks, and that the sea quarks do not contribute much to these processes (or that they are 
generated by higher orders). One can then show that quark-quark scattering can be simply 
embedded into the proton |9|, ||. The elastic amplitude for two-gluon exchange between 
two protons has been shown to be: 

A 2 = 8ta 2 s J d 2 Kd 2 k b J^l^ MK + k 6 ) - S 2 (K, k fe )] 2 (2.6) 

Si and £ 2 are two of the three form factors that can occur in the valence quark description 
of a proton. Each of these form factors correspond to the situation where 1, 2 or 3 quarks 
get hit by gluons. One can write their expression in terms of the proton wavefunction ip: 

£ 2 (KM) = J dM\^ j ,v j )\ 2 e ik ^ +ik ^ 
£ 3 (k a ,k b ,k c ) = JdM |^(/3 J ,r,)|V k - r '= +ik ^ +ikc - r - (2.7) 

where j3j is the fraction of longitudinal momentum (similar to Bjorken x), and r.,- is the 
transverse position of quark j. The natural integration measure dA4 is defined as: 

dM = [ I] ^A# (2) (E^ r ^(Eft - 1) (2-8) 

3=i,n q j j 

The first delta function defines the center of momentum of the hadron, whereas the second 
one enforces longitudinal-momentum conservation. Assuming that hadrons are made of 
valence quarks only, we normalise the wavefunction according to: 

£(0) = |dJW|V>(&,r,)| 2 =l (2.9) 

These expressions become useful once one realises that the form factor £\ also occurs in 
the elastic 7*p cross section and is none other than the Dirac elastic form factor in Eq. ( |2.3|) . 
Hence one of the form factors is determined. Furthermore, £2 and £3 are related to £\ by 
the following properties: 

£ 2 (k a ,0) = £x(k a ) 
£ 3 (k a ,k 6 ,0) = £ 2 (KM) (2.10) 

These properties ensure the IR finiteness of Eq. (|2.6|) . One can either calculate these form 
factors using a model for the proton wave function |§, or use a simple parametrisation, 
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which takes into account both the infrared and the symmetry properties of the form factors: 
one simply needs to make an ansatz for £ 3 and obtain £ 2 from Eq. ( p.lOj) . We take 

£ 3 (k a , k b , k c ) = E x (k 2 a + k 2 b + k 2 c - c{k a ■ k b + k a • k c + k b ■ k c )) (2.11) 

where c is an arbritrary number which can be shown to be of order one in the proton 
case |J. 

We have written the form factors as functions of the gluon momenta to manifestly show 
that Eqs. (10) hold as the gluon momenta go on shell. However the form factor comes from 
the proton wavefunction and can only depend on quark variables. Hence k a , k b , k c should 
be thought of as differences between initial- and final-quark transverse momenta. 

At this stage, one is in a situation to describe elastic scattering and hence, through the 
use of the optical theorem, total cross sections. The model contains two free parameters, 
as and c, which we can tune to reproduce the data. One can then look at the shape of the 
differential elastic cross section to check whether the model makes sense. Unfortunately, it 
is well-known R |9], [15[ that the shape of the cross section comes out wrong: its logarithmic 
slope at the origin is infinite, and its curvature is too big. 

One may be tempted to look into higher order corrections for a solution to this prob- 
lem. Indeed, as is not small, and the perturbative series contains terms of order as +1 log n s 
which clearly can become big at large s. These terms can be resummed using the BFKL 
techniques but such formalism cannot account for diffractive scattering at low momen- 
tum transfers, where unfortunately most of the cross section is concentrated. The most 
obvious problem is that the rise of the cross section predicted by leading-log- s perturbative 
QCD is entirely different from that which is observed in data, as its leading contribution to 
the hadronic amplitude goes like s 1+2S5as . Secondly, the slope of the pomeron trajectory 
a' = 0.25 GeV -2 introduces a scale of the order of \j\fa' ~ 2 GeV, which comes in the 
description of the t-dependence of the hadronic amplitude. No such scale is present in per- 
turbative QCD, hence the differential elastic cross section retains its wrong shape. Finally, 
the amplitude is not factorizable: this is in fact a consequence of the infrared finiteness 
of the answer. Quark-quark scattering via gluon exchange diverges for massless gluons. 
Nevertheless, hadron-hadron scattering is infrared finite, as the colour of the hadron gets 
averaged for very long wavelength gluons. Hence there is a contribution that comes from 
the diagrams where gluons are exchanged between different quarks in the hadron. These 
diagrams feel the hadronic wavefunction, and hence their contribution depends on the 
target. 

Even the lowest-order two-gluon exchange cross section fails to reproduce the factoriz- 
ability of the hadronic amplitude. This prompted Landshoff and Nachtmann |l7j to pos- 
tulate that gluons have an intrinsic propagation length l//xo smaller than typical hadronic 
sizes. The propagator hence becomes: 

D(k 2 ) = d(k 2 /»l)/»l (2.12) 

with lim x ^o xd(x) = 0. In that case, the diagrams in which the gluons couple to different 
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quarks are suppressed by a factor ~ (25 GeV 2 )/Zq with respect to those in which the 
gluons couple to the same quark line, and factorisation can be restored, for // ~ 1 GeV. 



This kind of picture has recently been questioned on theoretical grounds jL8), but it 



seems to remain the only one available which reproduces the observed properties of the 
pomeron. In the following, rather than worrying about the ultimate nature of the pomeron, 
we shall adopt the following pragmatic approach: we are after rates for the production of 
rare particles produced in a rapidity gap. This is very similar to single-diffractive processes 
and elastic processes, except that we need to insert a rare particle production vertex in 
the diagram. The best phenomenological description of such processes is two-gluon ex- 
change, multiplied by a pomeron s-dependence given by Eq. Q2.1J) . We shall first set up 
the lowest-order perturbative calculation of Higgs production via two-gluon exchange. We 
shall then evaluate the importance of the IR region by using a constituent-gluon propagator 
which matches to the perturbative one at high k 2 ||19|| . The model hence contains several 
parameters: as, c and a non perturbative scale /io- We shall tune these parameters to the 
highest-energy data available, reproducing the total, elastic and diffractive cross sections at 
the Tevatron. Hence our calculation of Higgs production at the LHC will be parameter- free. 

Note that this procedure also factors in the gap survival. Indeed, diffractive cross 
sections surely produce a gap, and the final-state multiple interactions which may spoil it 
are taken into account by our tuning of the parameters. As the kinematics of the process 
we are considering is rather similar to diffractive scattering, we expect that the choice of 
parameters which accounts for gaps in diffractive scattering will also lead to gaps in Higgs 
production. 



2.2. A QCD model for inelastic diffractive cross sections 

Before embarking into the afore- mentioned fit, we need to develop the formalism needed to 
describe inelastic diffractive cross sections, i.e. those in which a proton gets hit by a colour 
singlet, but nevertheless fragments. The main ingredient needed is a new form factor for 
the proton. 

Hence we want to describe the process p — > X instead of p — > p. The simplest way to 
do this is to square the amplitude, and to recognise that the cross section, including all 
the interference terms in the final state, is similar to the elastic amplitude resulting from 
4-gluon exchange, with gluons arranged in two singlets. 

To calculate the corresponding form factor, one has to notice that one-photon exchange 
leads to an E\ form factor. The two-gluon form factor can then be built by squaring the 
one-photon exchange diagram, and removing the intermediate proton. This can lead to an 
Si form factor, if the same quark is hit in the two interfering diagrams, or to an £ 2 form 
factor, if different quarks are hit. Similarly, an £ 2 convoluted with an E\ will give either 
an £ 2 , or an £ 3 . In the case of the calculation of the elastic form factor, there is an extra 
complication coming from the colour algebra. However, in the case at hand, the colour 
algebra is trivial, because we are exchanging colour-singlets, and hence the colour factor 
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Figure 2: One of the 81 interference terms contributing to the inelastic diffractive form factor 

is factorised between the two sides of the diagram. Therefore, after we recognise that the 
elastic form factors are given by 9 terms, three £\ and six £2/2, we can simply write down 
the 81 possible interferences, and decide by inspection which form factor will describe them. 

More formally, one first needs to assume that in diffractive scattering one goes through 
an intermediate state, call it Y, described by a wave-function ipy The only change in the 
preceeding formulae (|2.7| ) then occurs in the substitution 

rj)\ 2 — > ip*(j3j,rj)ipY(/3j,rj). This means that we are not going to describe interme- 
diate states which have a mass very different from that of the proton, as their production 
necessitates a change in the longitudinal momentum. Hence the cross section we are calcu- 
lating represents the bulk of diffractive scattering, but does not reproduce the high-mass 
tail. This tail presumably comes from higher-order diagrams, be they 3-pomeron vertices 
or gluon radiation |2(| |P2| . 

Now, in elastic scattering, the form factor describing p — > p is 3(£i(t) — £i{t))- This 
comes from 9x2 diagrams similar to those of Fig. [1], as each gluon can be attached to 3 
different quarks. The form factor describing p — > Y is the similar, with e.g. £\ given by 

£x(k) = J dMi(>;(P j ,r j )MPj,r j )e ik - e '' (2.13) 

Squaring the amplitude, we obtain the square of these form factors. Let us for instance 
concentrate on the £\ x £ 2 term shown in Fig. |2|. It corresponds to the expression: 

^(k+ko^K^K) = E/tn ^^(EtoWE/W) 

Y Y j=l,n q j j 

r^r^i, r i )e i(k+kl) ' ri 
x [n«iW (2) (EW(E B rl) 

l=l,n q I I 

r Y (B u R0W(A, R,)e iKl - Rl+K - R2 (2.14) 
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The completeness relation for the wavefunctions reads: 

/(A,r<) = e/[ n ^R# (2) (E^-wE^ - 1) 

Y j=l,n q j j 

MPh nW Y (B 3 , R;)/ Rj) (2.15) 

Hence formula ( J2.14 ) becomes: 

= £ 2 (k + k x + K x , K) (2.16) 



Y 



All the other 80 interference terms can be worked out similarly, and we get the form factor 
in the inclusive case from the square of the amplitudes. 

Noting that the four momenta entering the form factor have to sum to zero (as we 
are squaring an amplitude, the initial and the final states must be identical), we use the 
short-hand notation: £2(1) = ^(l) — I), £3(1? L) = £ 3 (1, L, — 1 — L). We then obtain the 
following: 

3 T (k,k 1 ,K 1 ,K)=3[l-£ 2 (k)-f 2 (k 1 )-f 2 (K 1 )-f 2 (K) 

+ 2£ 2 (k±k 1 )±~£ 2 (k±K)±~£ 2 (k±K 1 ) (2.17) 

- £ 3 (k, kx) - S 3 (K, K x ) + ±£ 3 (k!, Kx) + ^ 3 (k, K x ) + ±£ 3 (k, K) + ^ 3 (K, k x )] 

Note that we recover the same IR behaviour of the square form factor as we had for 
its elastic components: jF(k 1; k 2 , k 3 , k 4 ) — > for any k.f — > 0. This insures the infrared 
finiteness of the answer. 

The diffractive cross section then has the same form as the square of the elastic ampli- 
tude, except that the form factor 81 [Si — £ 2 ] 4 gets replaced by 9JF 2 . 



2.3. Best values of parameters 

We are now in a position to fix the parameters of the model. As we need to extrapolate 
to the LHC, we use Tevatron data. The fits we obtain are good to about 10% in the 
perturbative case. Slightly better results are obtained if we smooth the infrared region of 
the propagator, using a constituent-gluon propagator, as explained at the end of Section 
2.1. 



Table 1: Diffractive data 



experiment 


o tot (mb) 


a ei (mb) 


2a sd (mb) 


CDF 


80.03 ±2.24 


19.7 ±0.85 


9.46 ±0.44 


E710 [§ 


72.2 ±2.7 


16.6 ±0.7 


9.37 ±2.9 
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Figure 3: The Feynman diagrams contributing to H boson production in a rapidity gap. 

The E710 data are well reproduced by the following set of parameters: 
Table 2: Parameters reproducing the data of Table 1 



propagator 


a s 


c 


Ho (GeV) 


perturbative 


0.88 


0.61 




constituent [|19| 


1.52 


0.11 


1.9 



Note that in the case of the constituent gauge propagator, the running of the coupling 
is included in the gluon propagator, and that in this case the value of as is that at the 
renormalisation point /i . 

3. Heavy Particle Production 

The purpose of this work is to estimate the cross section for the production of a heavy 
particle H within the central region, hence the heavy particle must be produced within the 
gluon exchange, and not as a result of the proton fragmentation. One can describe the 
production mechanism via an effective gluon-gluon-H vertex: 

, _ _ on Wo 

W = d^ih ■ k 2 g» u - Kk v 2 )-± + {k\klg» v + kltth ■ k 2 - k^k\kl - KKkX)—^ (3.1) 

and the specialisation to a specific particle is done via the calculation of W\ and W 2 - 

Once this vertex is known, we must embed it into the two-gluon exchange graphs. One 
can couple the vertex to either gluon, and there are two two-gluon exchange graphs, hence 



10 



V 


1 <2i 

1 


i PI 


^ i 
— i » 


p 

/ 


> 


v 1 > 






^ » 






oh 1 o 

r— j— -g 










ot^2 1 ^2 to 














— T~ 

p 


92 


i P2 


!<?2 


» * 

1 

p 



Fi gure 4: The square of the Feynman diagram contributing to H production in a rapidity gap. The 
imaginary part of this diagrams gives the total cross section for the process in Fig. pj. The dash-dotted line 
shows the cuts used in obtaining the imaginary part. 



in principle we have to calculate four diagrams. However, as noted in Ref. |I(J, the sum of 
the four diagrams in Fig. ^| is equal to the s-channel discontinuity of the first one. Thus we 
can obtain the quark-level H production cross section by calculating the imaginary part of 
the diagram in Fig. f|. 

Heavy particle production forces the gluons to carry a non-vanishing longitudinal mo- 
mentum, hence the kinematics is not purely transverse anymore. We shall have to modify 
slightly the expressions for the form factors to take this into account, and these form factors 
will then allow us to go from the process qq — > Hqq to pp — > Hpp, or pp — > HX1X2. 



3.1. Kinematics 



The expression for the quark-level cross section is given by: 

da = ^\M\* (3.2) 

with the 2/9 the colour factor for two-gluon exchange. 

The contribution to the total cross section of the process depicted in Fig |] is arrived at 
by calculating the imaginary part of the diagram in Fig f| via cutting rules. In the notation 
of Fig. H the square of the imaginary part of the amplitude given by: 

\M\ 2 = ^1 ^d 4 KdV^ 2 5 + (^)5 + (p 2 2 )5 + (gD5 + (g 2 2 )5 + (g?)^(Q2) 

5+(H 2 - M 2 H )M (3.3) 

where Ai is the squared differential amplitude given by Feynman rules (see the next section) 
and 5+(k 2 ) = S(k 2 )9(k ). 
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Let us rewrite the momenta in terms of Sudakov variables: 



k = 


x y . 
~P+ ~P + k 




s s 


Pi = 


, Vi / , 

XlP + — p + Pi 




s 




%2 / 


P2 = 


—P + V2P + P2 



with k, Pi, p 2 transverse to p and p'. 

We assume that the off-shellness of the incoming quarks can be neglected, so that 



%1 2 ) 


= % + k 2 ) 




= 5(x-k 2 ) 


s(pI) 


= ^5(y 1 + ^) 


x 1 x 1 


s(pD 


= -5(x 2 + ^) 

V2 V2 



(3.4) 
(3.5) 
(3.6) 



p.p = 0,p.p' = s/2,p'.p' = (3.7) 
One can then solve for the 5 functions putting the quarks on-shell: 

(3.8) 
(3.9) 

(3.10) 
(3.11) 

Treating the 5 + {Q\)5 + {Q 2 ) in a similar manner, we arrive at: 

d 4 k d 4 K 5{ql)5{ql)5{Ql)5{Ql) = -^d 2 kd 2 K (3.12) 
d 4 Pl d% 5(pl)5(pl) = ]^^d 2 Pl d 2 p 2 (3.13) 

4 X\ y 2 



We now have to deal with the (^-function putting the H on-shell: 

5((p 1 +p 2 -p-pT-M 2 H )=5(2p 1 .p 2 + s(l-x 1 )(l-y 2 ) + ^ + rf + V 2 _ M 2 ) {3U) 

x±y 2 X! y 2 

We will eventually use it to eliminate X\ or y 2 . 
Hence the differential cross section becomes: 

da = 1 , a j ^ ^d 2 p ld 2 p 2 d 2 kd 2 K 

81 x 8s (2tt) 5 s 2 xi y 2 H H 

x 5(s(l-x 1 )(l-y 2 ) + (p 1 + p 2 ) 2 + ^- + ^-M 2 H ]M (3.15) 
V xi y 2 J 

The positivity of the energy of the on-shell lines dictates the integration limits: 

0<x 1 ,y 2 <l, \y + x\<s, |F + X|<s. (3.16) 
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All the dot products can then be expressed in terms of the transverse kinematics and 



of the longitudinal momentum fractions x\ and y 2 . The integration bounds (|3.16| ) do not 
guarantee the existence of a rapidity gap, but simply that of a physical process. 

Several conditions are necessary for the existence of a gap. First of all, the final state 
must not be too far off-shell. If the incoming proton has momentum P, then the outgoing 
X cluster has momentum P x = P — P + Pi- Its mass is given by (write P = p/xb = ^bP) 
M\ 



\y>i\ 2 (u)b — 1)/ X\. Similarly for the outgoing cluster X' with momentum Vx' 



V -p' + p 2 we have M x , = \p 2 \ 2 {uj' B ~ 1 )/V2- 

The rapidity gap Ay between the two outgoing quarks is given by: 



Ay = § In 



P -Pi 
p.pi 



I In 



P-P2 
P-P2 



hi 



IP1IIP2I 



and that between the two hadronic clusters is: 



Ay 



V'.Vx 



Mn 

2 V.Vx 



v.v x > 



h In 

2 P.Px' 



log 



2 xiy 2 (uj B + xi- l)(u' B + y 2 



(3.17) 



(3.18) 



Setting aside the question of gap survival (see section 4.4), a large rapidity gap will be 
present if the fractional longitudinal momentum loss of the quarks 1 — xi, 1 — y 2 , and the 
transverse momentum fraction, |pi||p2|/s, are small. The gluon propagators automatically 
provide us with small |p>i| |p»a|/ s- 

Whereas large gaps are produced for x\ and y 2 near one, heavy particle production 
requires that neither X\ nor y 2 be too close to unity: 



Ml 



(1 - £i)(l -y 2 )s 
Since for valence quarks, s ~ s/9, we arrive at the kinematic limit 

M H < (l-xOv^/3 



(3.19) 



(3.20) 



We shall take X\ > 0.7, since cross sections then become mainly diffractive [12]. Roughly 
speaking, particle production and detection in a rapidity gap is limited to masses less than 
an order of magnitude below the center of mass energy. Note that here x\ is defined at the 
quark level, and that this cut is similar to that used by Bialas and Landshoff. 



3.2. The s-channel discontinuity of the amplitude 

Let us now turn to the evaluation of \M.\ 2 (see Eq. ( |3.3|) ), which is what couples to the 
diffractive form factor. Because of the flat high-s behaviour of the cross section, and 
because the exchange is even under charge parity, the amplitude M. is equal to its s-channel 
discontinuity. 
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M. is made of two traces, T\ and T 2 , six gluon propagators, D(k 2 ), and two effective 
vertices from the Higgs loop, $f^ A and <fi± x . We thus have, (see Eq. (|3.3|) ): 

M = T 1 T 2 ^ K <p\ x D{kl)D(k 2 2 )D(k 2 )D(Kl)D(Kl)D(K 2 ) (3.21) 

The traces are given by: 

T\ = Tr(p.j 7^ gi.7 7^ p 1 .j 7^ Q1.7 7^) 

T 2 = Tr(p'.7 7^ g 2 -7 7a P2-7 7a <?2-7 7a/) (3.22) 
and the tensor structures in Eq. (|3.1|) associated with the ggH vertex are: 

tt X = {k 1 .k 2 g vX -k^)W 1 /mj I (3.23) 

Here, we have neglected the W2 term, as it will be smaller in the special instance of Higgs 
production (see section 4 and the Appendix). 

The number of terms arising from the traces, after we perform a Sudakov expansion of 
the momenta, is of the order of one million and so we will keep only the leading order in 
the quark center of mass energy s. The final answer then takes the form: 

M = s i AD(k 2 l )D(kl)D(k 2 )D(Kl)D(Kl)D(K 2 ) (3.24) 

Here D is the gluon propagator and A is a complicated function of the transverse momenta, 
which includes the triangle vertex. We do not quote the full expression for A as used in the 
program since that is not very illuminating. Instead we present the first terms in a series 
expansion in 1 — xi, 1 — y 2 , which as discussed in the previous section, are required to be 
small (< 0.3) to insure a large rapidity gap. 

A=(A + Ai(l - an) +A 2 (1- y 2 )) WWfA4 (3-25) 

with: 

Ao = 64 (k. Pl - k.p 2 - k 2 ) (K. Pl - K.p 2 - K 2 ) 

At = -64 [(k. Pl - k 2 ) (2 K. Pl - K.p 2 - 2 K 2 ) - (K. Pl - K 2 ) k.p 2 ] 

A 2 = 64 [(k.p 2 + k 2 ) (K. Pl - 2 K.p 2 - 2 K 2 ) + (K.p 2 + K 2 ) k.pj (3.26) 

3.3. The form factor 

We have discussed elastic and diffractive form factors when all gluon momenta are trans- 
verse. We are however interested in producing a Higgs particle in the rapidity gap. This 
requires that a significant fraction of the proton longitudinal momentum be carried by the 
gluons. Our previous expressions for the form factor needs to be adapted to this situation. 
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Figure 5: Valence quark description of proton where one (a) or two (b) quarks get hit by a gluon. 



Let us first consider the process pp —> ppH, similar to the elastic cross section in 
Eq. (|2.6|) . This involves two form factors corresponding to diagrams (|5]a) and (||b). One 
starts with a proton of momentum V and ends with an on-shell proton of momentum 

2 

V 1 = x{P - —V + Pl (3.27) 

XiS 

V' is the momentum of the proton with three momentum opposite to that of V. What 
matters for the quark recombination, and hence for the form factor, is the transverse 
kinematics in the V\ frame. The total momentum transfer is A = V\ — V — k + k\, with 

k = 0(l/s)V + 0(l/s)V' + k 

h = (xt - l)V + 0(l/s)V' + pi - k (3.28) 

We want to calculate the transverse parts of the gluon momenta k\ and k w.r.t. the V\ 
direction. The momentum of the first gluon k is purely transverse to order 1/s and does not 
directly enter into Higgs production. To order 1/s its transverse part remains unchanged. 
The momentum k\ has a non-negligible component in the V direction. With respect to the 
V\ = V + k + k\ direction it becomes 

h = (^— ^)CP + k + h) + ki + 0(l/s) (3.29) 

Xi 

where the vector k x is now transverse to V\. Solving this equation for k\ we get 

^ = (x 1 -l)P + 0(l/s)P' + ^^k + k 1 (3.30) 

X\ X\ 

which by comparing to Eq. (|3.28p gives: 



k 1= Pl^-^k=Pi-k (3.31) 

X\ X\ X\ 

Hence, the transverse momenta which must enter the form factor are k and ki. We need 
to worry a little further as the coefficients of the dot products can also be x\ dependent. 
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To check this, and to check how x\ enters into the diagram, we consider the simplest 
case: the form factor £ 1 in the case of elastic of pion scattering, way is to consider £ 1 in 
the case of pion exchange. 

First solve for the pion kinematics: let A = (x± — 1) V + 2/3a/ sP' + A, with V.V = 0, 
V'.V = and V.V = s/2. The on-shell condition for the final pion, (V + A).(V + A) = 0, 
gives: 




(3.32) 



and the total momentum transfer squared is: 

A 2 

t= — =an(ki + k) (3.33) 

X\ 

As we are calculating the imaginary part of the diagram, we have a few on-shell condi- 
tions for some of the quarks: one of the incoming quarks is on-shell, call its momentum r, 
with r = x V + 2[3 r /s V + r. The on-shell condition r.r = gives 




(3.34) 



The off-shellness of the other quark is: 

2 

m 2 = (V -r).CP -r) = — (3.35) 

x 

The momentum of the first gluon is given by k = 2a k /s V + 2f3 k /s V + k, and the last 
on-shell condition, (V — r + k).(V — r + k) = 0, gives: 

k 2 x — 2 k.r x + r 2 _ . , . . 

ft= 2,(,-l| + °'^» < 3 - 36 ' 

The second gluon is k\ — A — k, with 

h » ( Xl - 1) V + ~ 2 Xl/3k - Xl 1 V - k + A (3.37) 

X\ s 

Now, the dot product of the two quarks making up the proton is: 

A 2 = 2r.(p - r + A) = -2(p - r + A).(p - r + A) = -m 2 xi + 2r.A - x t (3.38) 

The form factor emerges as a convolution of wavef unctions. These can only depend on 
the dot products of the 2 quarks making the proton, hence on A± and A 2 . In the case where 
Xi — 1, one gets: 




(3.39) 
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which clearly is a function of A only. In the more general case, one gets: 



f i)*C—)ii(^— - sfx~/x~ x Afd 2 Ydx (3.40) 

So, for t\ » m 2 , one gets F\{ti). However, this is not true in general, and the form factor 
clearly has an extra x\ dependence: A = does not lead to a form factor equal to 1. One 
can in fact make a rough guess: the effective argument is t = t+ < — > (1 — Xi) and 
< ™i >~ 9 x 0.3 2 ss 1 GeV 2 . 

X 

Hence we see that the ansatz is as follows: replace the transverse vector by ki, and 
multiply the overall argument by X\ . Note that this guarantees the infrared finiteness of the 
answer when k 2 or k\ = k 2 /xi — > 0. Clearly, this procedure straightforwardly generalises 
to the diffractive case pp — > XiX^H. 



4. Higgs production 

Apart from the form factors W\ and W%, the preceeding formalism can be applied to the 
production of any rare particle within a rapidity gap. We shall now specialise the calculation 
to the production of the minimal standard model Higgs boson. We first calculate the 
effective coupling resulting from the triangle graphs. 



4.1. Higgs-gluon-gluon effective vertex 

By calculating the diagrams in Fig. || we arrive at the following effective Higgs-gluon-gluon 
vertex: 



S^iki ■ k 2 g^ - KK)—^- + (klklg^ + m&kx ■ k 2 - k^k'fkj - k^k v 2 k\)—^ (4.1) 



where 



and 



Wi= W2G A ^ BNl (4 . 2) 
3tt 3tt 



Xi = 3 i dx f 1 X dy 1 ^ Xy (4.3) 



o Jo 



N 2 = ^f/'dx/" 1 ^ 1 - 2 "- 2 ^ 4 ^ (4.4) 
k\ ■ k 2 Jo Jo D 

D = \-{2k 1 -k 2 /m 2 t )xy + {kl/m 2 t )(y 2 ~y) + (kl/m 2 t ){x 2 -x) (4.5) 
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— - + 





Figure 6: The Higgs-gluon-gluon vertex. All quarks run around in the loop, but since the Higgs-quark- 
quark coupling is proportional to its mass, the top quark contribution dominates. 



Because of the kinematics of Higgs production in a rapidity gap, k\ , k\ << 2k\-k 2 
and we can approximate the integral expression by setting kf, k\ equal to zero: 



m 



N° 



N° 
D° 



3 1 dx r x d, }~ Axy 

6 Jo Jo V D° 

a f 1 a t l ' x a 1 ~ 2x - 2 y + 4x y 

b ax ay 

Jo Jo D° 

m 2 H / m 2 )xy — it 



1 



(4.6) 

(4.7) 
(4.8) 



Note that with k\ and k\ set equal to zero it is the ie prescription coming from the top 
quark propagators in the loop that defines the contour. 

Since the off-shellness of the gluons is small compared to the Higgs mass, and since the 
tensor structure (|4.1|) will always contract k\ with the upper side (//) of the diagram and 
&2 with the lower side (u), the W2 term in Eq. ( |4.1|) can be ignored, unless of course \N$\ 2 
turns out to be abnormally large compared to |iV°| 2 . As we explain in more detail in the 
Appendix, \N% | 2 is always at least 30% smaller than \N°\ 2 for Higgs masses of 1 TeV or 
less. Thus we ignore W2 and the above effective vertex can be derived from the following 
momentum space Lagrangian 



£ eff = A a tl A av H{kx ■ k 2 g^ - Kk v 2 )^ 



in 



(4.9) 



H 



The term above, while not gauge invariant, forms part of the following gauge invariant 
term. 

1 , ; „ ,, (1.10) 



eff 



2m 2 ^ ^ 



As a check we use our effective vertex to calculate the decay for H — > gg and we obtain 
the well-known result 

IH^I 2 



m - gg) = 2 , AixmH 
where the last 2! is an identical final state particle phase space symmetry factor. [ 



(4.11) 



3 Our result for the effective vertex differs from the result quoted by Bialas and Landshoff by a factor 2, 
presumably the above symmetry factor. 
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The approximate form of N\ in Eq. ( [4.8| ) can be evaluated in closed form. Defining 
a = m 2 H /m 2 we have 



N? = 3 f 1 dx f 1 " dy 1 AXV 

Jo Jo 1 — axy — le 

= - 1 - (1 - 4/a) ^arctan[y / a/(4 - a)] 
= - 1 + (1 - 4/a) ^arctanh[ya/( a - 4)]^ 



if a 



m 



mi 



f < 4 (4.12) 



if a = 



m 



H 



ml 



> 4 



where arctanh(x) = In \ (x + l)/(x — 1)| — z7r/2. At a = 4, iV° = 3/2, and iV° is continuous 
for the entire range of a. Note the change in sign in front of the second term between the 
two expressions. 



4.2. The Higgs production cross section 

From Eq. (|3.2j) we get: 



s J 8 (2^)^ j x±y 2 

&F(Jfe, ki, K x , K) 3f(k, k 2 , K 2 , K) (4.13) 

with k 2 ^ (k + Pl ) 2 + (1 - xi)(k 2 - p 2 ), A; 2 « (k + p 2 ) 2 + (1 - y 2 )(k 2 - p 2 ), A; 2 « k 2 , and 
similar expressions for K 2 , Kf, fC 2 , and the Higgs momentum 

2 2 

H 2 = ( Pl +p 2 -p- p') 2 = s(l- Xl + ^)(1 - y 2 + + (pi + P2 ) 2 (4.14) 

sy 2 sxi 

The form factors 3J r (k, fci, i^i, K) 3jF(fc, fc 2 , if 2 , ^0 correspond to calculating the inelas- 
tic cross section pp — > HX\X 2 . To calculate pp — > ppif these factors are replaced by an 
expression written in terms of the Dirac elastic form factor Fi (t) given in Eq. (|2.3| ) 

(3Fi(ti)Fi(t 2 )) 2 where t x = p 2 /xi and t 2 = p 2 /y 2 (4.15) 

In order to make explicit the flat s behaviour we first use Eq. ( |3.24|) to substitute for M. 

o = / ^^d 2 Pl d 2 p 2 d 2 M 2 K5(H 2 - M 2 H )AW X W[* /m A H (4.16) 

18 (ztt) J x\y 2 

F(k, h, K 1} K)F(k, k 2 , K 2 , K)D{kl)D{kl)D{k 2 )D{Kl)D{Kl)D{K 2 ) 

J\4 brings in a factor of s 4 and we still need a factor of 1/s so we solve the Higgs momentum 
5 function for x\ or y 2 . (Note that scaling pi and p 2 by yfs produces the wrong power of 
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s.) In order to facilitate presenting the solution, we linearise the kinematics in terms of 
1 — y 2 . Hence 



2 ^2 ^2 ^2 



H 2 » 5 ((i_ Xl ) + P2 + (i_ y )2)^)((l- 2/2 ) + ^ + (l-x 1 )^) + ( Pl + p 2 ) 2 

s s s s 

« S((l-x 1 ) + ^)((l-y 2 ) + ^) + ( P i + P2) 2 (4.17) 



Thus solving the 5 function we get 

i-^d-ri ^r'v^ - < 4 ' i8) 

This brings a 1/s/ ((1 — Xi) + p|/s) from the 5 function upon performing the eh/2 integral 
and the answer behaves like s°. 

There is an overall arbritrary angle in the integral with respect to which we measure 
all angles. We choose to measure our angles from k. Thus d 2 h = d8dh 2 /2 = ndk 2 . Thus 
we arrive at: 



a 



<4 



J d%\ d\a d 2 pi d 2 p 2 d 2 ~K 



x FfckuKuK) F{k,k 2 ,K 2 ,K) (4.19) 
x D{k{) D(k 2 ) D{k 2 ) D{Kl) D{Kj) D(K 2 ) 

As explained in section 2, we need to multiply the above expression by a Regge factor, 

rii=i,2 (f:) ) with Si = (pi + H) 2 , and U = (p — pi) 2 . This is our final expression for 
the cross section. Note that although the Regge factor is very important when one fits the 
purely soft cross sections (as in section 2.3), it plays very little role here, because M\, and 
hence Si are sizeable fractions of s. 



4.3. The Bialas-Landshoff Limit 

Bialas and Landshoff (BL) [l(J work out the leading term in (1 — xi), (1 — y 2 ) at pi = 
P2 = 0, and then reintroduce the pi, p 2 integration after reggeisation. In this limit, we 
get: k\ = Xik 2 , k\ = y 2 \s?, A = (8 |k 2 |Wi) 2 s 4 . To compare more easily with BL, we do not 
eliminate the on-shell Higgs condition, and write: 

, 4 a 4 dx 1 dy 2 , 2 , 2 P//1 N/1 . Mfj 

do = g1 - „„, * -d 2 Pl d 2 p 2 S {{1-xt I-2/2 --^ 

81 x 8s (znys 6 Xi y 2 s 



X 



[/ rf 2 k8 \k 2 \D(k 2 )D(x 1 k 2 )D(y 2 k 2 )W 1 ] 2 s 4 
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a s dx l dy 2 2 2 „ . Af£. 

51^ X\ 1/2 S 



X 



[/ d 2 k \k 2 \D(k 2 )D(x 1 k 2 )D(y 2 k 2 )W 1 } 2 (4.20) 



which is identical to formula (4.2) in Ref. [HJ at e = a' = A = 0, except for the factor of 2 
which come from symmetrizing two identical particles as discussed in section 4.1. 

Note that this is not the only simplification in their work, as they also assume that an 
exponentially falling propagator multiplying F\ is a good representation of the product of 
the form factors by the gluon propagators. Furthermore, their estimate of the diffractive 
cross section comes from the replacement of Fi by 1, and their integrals then converge only 
because of the Regge slope. We shall now see that, surprisingly, the improvements we have 
introduced - higher orders in (1 — xi), full form factors, propagators - do not significantly 
modify the final answer. 



4.4. Results 

Eq. (|4.20| ) is an eight-dimensional integral which we evaluate with the Monte-Carlo program 



VEGAS p2|, both for the elastic cross section, pp — > ppH, (corresponding to using the 



form factors in Eqs. (|2.7|) and ( |2.11|) ), and for the diffractive cross section, pp — > X H, 



(corresponding to the form factors described in Eq. ( [2.17|) ). 

We present our results for the Tevatron (y/s = 1.8 TeV), and the LHC (\/s = 10 TeV 
and y/s = 14 TeV), in Figs. f7LpL|9|. The bands represent the effect of switching from a purely 
perturbative propagator to a renormalisation group improved propagator suppressed in the 
infrared region. We see that the extrapolation from the purely diffractive cross sections of 
section 2 to Higgs production seems to bring in relatively modest uncertainties of about 
a factor 2. The upper bands are for diffractive scattering, the lower bands for elastic 
scattering. 

It is striking to note that our calculation leads to essentially the same conclusion as 
the Bialas-Landshoff estimate. This can be understood for the following reasons. One 
phenomenologically parametrises the model so as to correctly reproduce the soft elastic 
and diffractive cross sections. In the Higgs production diagrams, the off-shellness of the 
gluons is rather small < 1 GeV 2 , and comparable to that in the soft diffractive processes 
which we fit to. The deep infrared region is cut off by the form factors. Thus the behaviour 
of the gluon propagator in the UV region and the deep IR region does not matter much, 
and the uncertainties arising from the extrapolation to Higgs production are small. 

The diffractive cross section (upper band) is on average about a factor 5 higher than 
the elastic one (lower band). We call attention to the fact that it is really a factor 10 
higher in the purely perturbative case, and a factor 3 higher if we use the constituent-gluon 
propagator. In other words the perturbative calculation is the upper limit of the diffractive 
band and the lower limit of the elastic band. Given the smallness of the gluon off-shellness, 
we consider the perturbative result less reliable. 



21 



.6 



.4 



b 



.2 



.0 




50 100 150 

M H (GeV) 



Figure 7: The Higgs production cross section within a gap at 1.8 Tev 
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Figure 8: The Higgs production cross section within a gap at 10 TeV 
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Figure 9: The Higgs production cross section within a gap at 14 TeV 



24 



We see that the cross section at the Tevatron is rather small, of the order of a fraction 
of a picobarn, for light Higgs masses. It indeed gets rapidly cut by the condition (|3.20|) 
which precludes production of particles of mass bigger than 200 GeV. If the Higgs is below 
100 GeV, it will be observable within a rapidity gap at the Tevatron. A handful of them 
may already have been produced within the 60 pb -1 available today. 

Fig. [10] helps us understand the shape of the Higgs production curves, Figs. |7]||[|. 
Fig. [lC)| a shows the vertex as a function of the Higgs mass. It rises with the Higgs mass, until 
one reaches the threshold for top quark production (after which the vertex is dominated 
by its imaginary part). Fig. |T0| p shows the x± distribution for a 100 GeV (solid line) and 
a 550 GeV Higgs (dashed line), in the diffractive case, for y/s = 10 TeV. Although the x\ 
distribution is higher for the higher mass (because the vertex is bigger), it gets cut off at 
x\ ~ 0.9, after which one does not have enough partonic energy to produce a 550 GeV 
Higgs. Hence the low-mass Higgs cross section is higher because it is dominated by large 
X\. This also shows that the dependence of the cross section on the cut-off X\ > 0.7 (used 
to insure a rapidity gap) is very weak at low Higgs masses, and gets progressively larger 
as the cross section decreases: the main cause of its large fall-off at large Higgs masses is 
indeed the X\ cut. Putting these considerations together we can understand the dip in the 
Higgs production curves before the rise peaking at about twice the top mass. Fig. |TT| gives 
an estimate of the size of the rapidity gap. We show the rapidity distribution of the Higgs 
(thick solid line) and the partonic cluster (thin solid line). As the latter will fragment into 
hadrons, we envisage a worst case scenario, where the cluster decays into a forward proton 
and a backward pion. The dashed curves give the pion rapidity distributions. A gap is 
clearly present for yfs = 10 TeV, and the Higgs peak becomes sharper as the mass of the 
Higgs increases. At the Tevatron, we see that the gap is reduced to a few units of rapidity, 
but remains present. 



5. Conclusion and further studies 

We have considered the lowest order QCD diagrams that give rise to Higgs production 
within a rapidity gap. They correspond to single pomeron exchange. We have found that 
the cross section is stable with respect to the details of the infrared region, and that it is 
a substantial fraction of the total Higgs production cross section this way at the Tevatron, 
whereas the LHC will comfortably produce 1 pb of 600 GeV Higgses within a gap. 

There are three main factors effecting the accuracy of our predictions: the background, 
the survival of the gap, and higher order corrections. The first two effects will reduce our 
estimate, the last will increase it. 

Although no conventional background is expected, there will be direct heavy quark 
production in the rapidity gap through the same color singlet exchange diagrams. Some 
work along this line |23| indicates that this background will be tiny for Higgs masses greater 
than 200 GeV. For Higgs masses less than 100 GeV the cc and bb background can be a 
problem. 
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Figure 10: (a) The Higgs production vertex as defined by Eq. ( 4.12 ). The thin solid line shows the real 
part, the dashed line the imaginary part, and the thick solid line the magnitude, (b) The X\ distribution 
for a 100 GeV Higgs (solid line) and a 550 GeV Higgs (dashed line), for \/~S = 10 TeV, in the diffractive 
case, using a constituent-gluon propagator. 
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Figure 11: The rapidity distribution of the Higgs (thick solid line), of the partonic cluster (thin solid 
line), and decay pions (dashed line) at (a) the LHC with ^/s = 10 TeV and a 100 GeV Higgs (b) the LHC 
with y/s — 10 TeV and a 550 GeV Higgs (c) the Tevatron with a 100 GeV Higgs. 
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Another issue to consider is the survival of the rapidity gap. Although the lowest order 
diagram produces a gap, further final-state interactions between the out-going partonic 
clusters might destroy it. Typical estimates for the survival of the gap give a number 
between 1 and 10% ||24|| . However, we want to point out that such an effect is implicitly 
included in our calculation: we fit the soft diffractive cross section, which already have to 
include such a survival factor, and as the kinematics of the event is not very different from 
the soft one, we expect those corrections to be small. 

Finally, we have set up the calculation in such a way that higher orders can be calculated 
la BFKL [[RJ. In the elastic case these corrections dramatically increase the cross section 
(e.g. a factor of y^/M^). We expect a similar effect in our calculation. 



6. Appendix 



In this Appendix we will justify neglecting the W 2 term in the effective vertex given in 
Eq. (}4.1| ) when calculating Higgs boson production in a rapidity gap. Given the cuts we've 
imposed to insure the survival the gap, the approximate form of N 2 given in Eq. ( |4.8j ) can 
be used. Defining a = m 2 H / m\ we have 

6 / dx 

Jo Jo 1 — axy — te 

-12 " 



N° 



i-z 1 - 2x - 2y + Axy 
dy — 



5 - (1 + A/a) arctan[Ja/(4 - a)] 



a/(4 - a) 



arctan[-i/a/ (4 — a)] 



if a = 



ni 



ii 



mi 



< A 



-12 



5 + ( 1 + A/ a) ^arctanh y a /(a — A) 



a/(a - 4) 



( arctanh-Ja/ (a — 4) 



if a 



m 



H 



mi 



> A 



(6.1) 



(6.2) 



(6.3) 



4, N° 



•15 + 37r 2 /2. It is the ie prescription that tells us how to analytically 



At a 

continue from a < A to a > A by replacing 



arctan 



A - a 



A - a 



i In 



a/(a - 4) + 1 
a/(a - 4) - 1 



i [ arctanh-/ I (6.4) 

V ya-A) 



We are interested in showing that the N$ is never much greater than the term. The 
analytic expression for is given in Eq. ( |4. 1 2| ) . We plot the ratio \N% \ 2 /\ N°\ 2 in Fig. |12 
and see that it is always less than 0.30. 
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Figure 12: The ratio \N% | 2 /|iV°| 2 as a function of a = m 2 H /m% between a = and a = 33. With a 
top mass of Tflf = 175 GeV this corresponds to a Higgs mass < rriff < 1 TeV. 

Acknowledgements 

JRC thanks A. Bialas, S. Brodsky and D. Soper for their comments and suggestions. 



References 

[1] HI Collaboration (T. Ahmed et al.) Nucl.Phys. B435 (1995); 

ZEUS Collaboration (M. Derrick et al.), Phys. Lett. B346 (1995) 399. 

[2] J.D. Bjorken, Phys. Rev. D47 (1992) 101. 

[3] P. Giromini, in the Proceedings of the V th Blois Workshop, Brown University, 8-12 
June 1993 (World Scientific: Singapore, 1994), p. 30. 

[4] R. Rubinstein, in the Proceedings of the V th Blois Workshop, Brown University, 8-12 
June 1993 (World Scientific: Singapore, 1994), p. 23. 

[5] G. Ingelman and PE. Schlein, Phys. Lett. 152B (1985) 256. 

[6] R. Bonino et al, Phys.Lett. 211B (1988) 239. 

[7] K. Goulianos, preprint RU-95-E-26 (May 1995), e-Print Archive: |hep-ph/950531C| . 
[8] J.F. Gunion and D. Soper, Phys. Rev. D 15 (1977) 2617 
[9] J.R. Cudell and B.U. Nguyen, Nucl.Phys. B420 (1994) 669. 



29 



A. Bialas and P.V. Landshoff, Phys. M.B256 (1991) 540. 

A. Donnachie and P.V. Landshoff, Nucl. Phys. B244 (1984) 322; B267 (1985) 690; 
Phys. Lett. B296 (1992) 227. 

K. Goulianos, Phys. Rep. 101 (1983) 169. 

A. Donnachie and P.V.Landshoff Phys. Lett. B296 (1992) 227. 

F.E. Low, Phys. Rev. D12 (1975) 163; S. Nussinov, Phys. Rev. Lett. 34 (1976) 1286 

E.M. Levin, M.G. Ryskin, Yad. Fiz. 34 (1981) 1114; 

L.V. Gribov, E.M. Levin, M.G. Ryskin Phys. Rep. 100 (1983) 1. 

E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45 (1977) 199 
Y.Y. Balitskii and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822 

P.V. Landshoff and O. Nachtmann, Z. Phys. C35 (1987) 405. 

K. Buttner and M.R. Pennington, preprint DTP-95-54, (June 1995), e-Print Archive: 
|hep-ph/ 9506449 . 



J. R. Cudell and D.A. Ross, Nucl. Phys. B359 (1991) 247 

J.R. Cudell, Proceedings of the 4 th Blois Workshop on Elastic & Diffractive Scattering, 
La Biodola, Italy (1991). 

M.G. Ryskin, Sov. J. Nucl. Phys. 52 (1990) 529, Nuclear Physics B (Proc.Suppl.) 18C 
(1991) 162. 

J.F. Gunion, H.E. Haber, G.L. Kane and S. Dawson, The Higgs Hunter's guide, 
preprint SCIPP-89/13, June 1989, (Addison- Wesley, 1991); Errata, preprint SCIPP- 
92-58, Dec 1992, e-Print Archive: |hep-ph/9302272l . 

G.P Lepage, preprint CLNS-80/447 (March 1980). 

A. Bialas and W. Szeremeta, Phys.Lett. B296 (1992) 191; W. Szeremeta, Acta Phys. 
Pol. B24 (1993) 1159. 

[24] R.S. Fletcher, Phys. Lett. B320 (1994) 373-376; 

J.D. Bjorken (SLAC). SLAC-PUB-95-6949, Jul 1995; 

E. Gotsman, E.M. Levin, U. Maor, Phys. Lett. B353 (1995) 526. 



30 



